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Abstract. A topological flux function is introduced to quantify the topology of magnetic 
braids: non-zero line-tied magnetic fields whose field lines all connect between two boundaries. 
This scalar function is an ideal invariant defined on a cross-section of the magnetic field, 
whose integral over the cross-section yields the relative magnetic helicity. Recognising that the 
topological flux function is an action in the Hamiltonian formulation of the field line equations, 
a simple formula for its differential is obtained. We use this to prove that the topological flux 
function uniquely characterises the field line mapping and hence the magnetic topology. A 
simple example is presented. 



1. Introduction 

In this paper we present a complete topological invariant for so-called magnetic braids: magnetic 
fields in a flux tube with monotone guide field. Such magnetic fields arise, for example in coronal 



loops in the Sun's atmosphere (Reale, 2010), or in toroidal fusion devices (Morrison, 2000). In 



both cases, the magnetic field is embedded in a highly-conducting plasma so that the magnetic 
topology — or linking and connectivity of magnetic field lines — is approximately preserved. As 
such, topological invariants typically play a significant role in the dynamics of the plasma and its 



magnetic field (|Woltjer[|1958[ |Taylor[|1986[ |Brown et a/.[|1999[ |Yeates et a/.[|2010[ |Candelaresi 
&; Brandenburg, 2011). 

Figure [1| shows our simply connected "flux tube" F, with a set (p, i\),z) of orthogonal 
curvilinear coordinates satisfying X7z • (Vp x V^) ^ 0. The magnetic field B is assumed to 
satisfy B z > everywhere in V, so that Vz aligns with the "axis" of the flux tube. We take ip 
to be an angular coordinate and p to be a radial coordinate in each plane of constant z. On the 
side boundaries of the tube, we require B • n = 0, so that all magnetic field lines stretch from 
one end Dq of the flux tube to the other end D\. We can then define the field line mapping 
F : Dq — >> D\ by F(xq) = /(xo;l), where /(x&z) denotes the magnetic field line rooted at 
x <E D . 

Two magnetic braids B, B are topologically equivalent if and only if one can be reached from 
the other by an ideal evolution 

— = V x (v x B) (1) 

with v\dv — throughout. We shall assume that their respective field line mappings match on 
the side boundary, F\q^ = F\qd , with the same winding number. In that case, B and B are 





Figure 1. Coordinates and notation for a 
magnetic braid B in a flux tube V, with 
field line mapping F : D$ —> D\. 



Figure 2. The physical meaning of Lemma 
[2] may be understood by integrating on the 
lower boundary Do along an open curve 7 
or a closed curve a. 



topologically equivalent if and only if F = F. Notice that F has two components. Our main 
result is to prove that a single scalar function, which we call a "topological flux function" is both 
necessary and sufficient to determine the topology. This we state in the following theorem. 

Theorem 1. Let B ; B be two magnetic braids on V , with respective field line mappings F, F 
that agree on 8Dq, with the same winding number. Let A, A be their respectivejtopological flux 
functions with the same reference field A ref such that A T ^ = 0. Then A — A if and only if 

F = F. 



It is well-known ( |Berger fc Field[ |1984[ |Brown et a/. , 1999) that having the same total 
magnetic helicity H r (which will be defined below) is a necessary condition for topological 
equivalence but not a sufficient one (two magnetic braids can have the same H r but different 
F). Indeed there are an infinity of other ideal invariants. For example, if XQi y Q (z) denotes the 
orientation of the line at height z between a pair of field lines f{x^z) and f(yo;z) : then their 
pairwise linking number 



-xo,xi 



1 

2^ 



de x ,yo( z ) 



dz 



is an invariant (Berger, 1986; Berger &; Prior, 2006[). So are analogous higher-order measures of 



the linking between triplets, 4-tuplets, etc., of field lines. Theorem [T] asserts that all of these ideal 
invariants are contained in A because, mathematically speaking, it is a complete invariant. It 



is the most economical description of field topology equivalence classes that is sought by Berger 



(1986). In future, the topological flux function should provide a useful tool not just to recognise 
but also to quantify changes in magnetic topology during non-ideal evolutions. For an initial 



application to measuring magnetic reconnection, see Yeates & Hornig (2011). 

For the alternative geometry of a half-space, Berger (1988) introduces A as the helicity of 



an infinitesimal flux tube around a single field line. Indeed, Taylor (1986) points out that each 
closed field line yields such an invariant in a perfectly-conducting plasma. This paper aims to 
interpret these constraints in a different light. We show that Theore m [l] is a consequenc e of A 
being the action the Hamiltonian system of the magnetic field lines ( |Cary fc Littlejohn , 1983). 
The Hamiltonian theory suggests an elegant expression for A in terms of differential forms that 
we exploit. 



A less general form of Theorem [I] has recently appeared in Yeates &; Hornig (2013). Here we 
present the results for a more general shape of flux tube with more general boundary conditions, 
and elaborate further on the Hamiltonian theory. We begin in Section [2] with the basic definition 
of A, before discussing the Hamiltonian interpretation in Section |3j The proof of Theorem [I] is 
given in Section [4] and two explicit examples in Section [5j 



2. Topological Flux Function 

We define the topological flux function A : Dq —> R as the line integral of A along the magnetic 
field line rooted at xq G Dq: 



rF(x ) 



A(x ) = / A(f(x ]z)) • dl. 

Jxq 



(3) 



As it stands, this definition is gauge dependent. Under a gauge transformation A — >► A+ V%, it is 
easy to see that A —> A + F*x~Xi where the pull-back notation means If*x)( x o) = x(^X x o))- 
We impose the gauge condition 



n x A 



dv 



n 



Aref I 
\dVi 



(4) 



where A ref is the vector potential of a reference field B ref that matches B • n on dV but is 
otherwise arbitrary. This particular choice of gauge condition is motivated by the relative 
magnetic helicity (Section |2.1[ ). Although this still leaves some freedom in x, & n d consequently 
in A, it is enough to make A an ideal invariant. To see this, assume an ideal evolution 



dA 
~dt 



vxBlVi (5) 

The boundary conditions of v|^y = and fixed n x A\gy mean that n x V$|ay = 0. Then, 
using the formula for the rate of change of a line integral over a moving domain (e.g., Frankel 



1997), we find 



dA 
d^ 



d_ r F 

dtJ X0 



F(x ) 



A dl 



/ f ^- - v x V x A + V(v • A) ) • dl 

F(x ) 

V(* + v • A) • dl = 0. 

XO 



(6) 

(7) 
(8) 



2.1. Relation to Magnetic Helicity 



Berger (1988) defines A as a "field line helicity", or limiting magnetic helicity in an infinitesimal 



flux tube around a single magnetic field line. Since our domain is magnetically open, we define 



a gauge invariant relative helicity ( Berger fej'ield , 1984; Finn & Antonsen, 1985) relative to the 
reference field B ref = V x A ref . Assuming Q, the relative helicity may be written 



(9) 



H r = [ (A + A ref ) • (B - B ref ) d 3 x = [ A • B d 3 x - [ A ref • B ref d 3 x. 
Jv Jv Jv 



The last term is a constant (call it i7 ref ) which depends only on the choice of reference field. For 
the first term, let xq be the footpoint of the field line through x = (p, -0, z), so that x = f{x$\ z). 



Changing coordinates gives 

H r - H Tei = f A(x) • B(x) d 3 x, 
Jv 

= f A(f( X0 ;z))-B(f( X0 ;z)) f z ^\ 
Jv ' B z {f{xo\z)) 

= / A(xo)B z (xo)d 2 x . 
Jd q 

Thus we see that A represents the helicity per field line. 



d 2 xodz, 



(10) 

(11) 

(12) 



3. Hamiltonian Interpretation 

We will show in this section that the differential of A has the following succinct formula in terms 



of differential forms (Frankel, 1997), where h p = ||<9(r)/dp||, h^ = \\d(r)/dt(j\\ are the coordinate 
scale factors. 

Lemma 2. Let B be a magnetic braid with field line mapping F and topological flux function 
A. Then 

dA = F*a — a, 

where a — A r ^h p dp + A^h^ dijj is the 1-form associated to A ref on Dq and (perhaps differently) 
on D\. 

This formula will be used to prove Theorem [l] in Section [4j It may be understood physically 
by integrating along a curve 7 G Do, so that 



f dA= f F*a- fa. 

J 'V i/'V i/'V 



In vector notation, equation (13) reads 



A{yo) ~ Axo) = I A • dl - / A 
J Hi) h 



dl, 



(13) 



(14) 



where xo,yo ^ ^o are the start and end points of 7 (see Figure [2]). If y$ ^ xo then 7 is an 
open curve, and (14) expresses the fact that the field lines rooted in 7 form a flux surface with 



§ A • dl = 0. This is the vertical shaded surface in Figure On the other hand, if y$ = xo 
then the curve is closed (the curve a in Figure [2]), in which case (14) expresses conservation of 
(vertical) magnetic flux in the corresponding flux tube. 



3.1. A as an Action 

It is insightful to derive Lemma [2] by interpreting A as an action in a Hamiltonian system. 
Suppose we parametrise a magnetic field line in terms of its length I as x(/), and think of the 
topological flux function as a functional, 



f F 

A(x ) = I 

JXQ 



F (xo) r l x 

A(x).-d<. 



(15) 



Then Cary &; Littlejohn ( 1983| ) point out that extremising A for given A and fixed end-points 
gives the path of the field line x(Z) through the domain. To see this, note that the Euler-Lagrange 
equations are 

8L A/ 8L ^ 0. (16) 



dxj dl \d(dxj/dl) 



In our case, the Lagrangian is L — A • (dx/d/), so 



dAi dxi dxi dAj 
dxj dl dl dxi 



= 0, 



which is equivalent to 



(Vx A) 



dx 
dl 



0. 



(17) 



(18) 



Hence dx/d/ is everywhere parallel to B, and x(/) is a magnetic field line. 

In fact, it is well known that the equations of the magnetic field lines are a Hamiltonian 
system. To show this, we follow Cary fc Littlejohn] (1983) and change the gauge of A so 
as to write the action explicitly in the canonical form of a Hamiltonian system. A gauge 
transformation A — > A + Vx changes the integrand (the Lagrangian) to L + dx/dl and also 
changes the integral (the action A), but leaves the Euler-Lagrange equations unchanged. If we 
set A p — everywhere in space, then A • dx = A^h^ dif; + A z h z dz. Making the identifications 
p «-»• A^(p, ijj, z)h^, q *-> ijj,t «-»• z,H *-> —A z (p, ij;, z)h z , our action becomes 



A 



/ V 

Jxo \ 



p— - H(p, q, t) J dt 



(19) 



which is a 1 degree-of- freedom Hamiltonian system in canonical form. The generalised coordinate 
is ^, the generalised momentum is A^/i^, and the Hamiltonian is —A z h z . Time is the z-direction, 
so our Hamiltonian is, in general, time dependent. 

The gauge where A p — everywhere may be called a canonical gauge since it renders the 
Hamiltonian system in canonical form. For Theorem IT] it is sufficient that n x A|^y be in 
canonical gauge, so in practice we restrict the gauge of A ref to be canonical, i.e. A 1 ^ = 0. The 
gauge A p = is equivalent to writing B in the form B = V(A^h^) x V^ + V(A z h z ) x Vz 



(Boozer, 1983; Yoshida, 1994). 



3.2. Proof of Lemma^ 

Given that A is the action in a Hamiltonian system, we can use known properties of Hamiltonian 

systems. Let f z be the "height-z" mapping of our field line system, defined by f z (xo) = /(#o; z). 

This mapping must preserve magnetic flux. In terms of differential forms, f*uo — uj — 0, where 

uj — B z h p h^ dpAd^. Every Hamiltonian system preserves such a symplectic form uj. But notice 

that 

U = (§p( A ^ ~ ^( A A)) dp Ad^ = da (20) 

where a = A^h^ d^ + A p h p dp. The primitive a is called the canonical/Liouville 1-form. 

Assume now that A is in canonical gauge as above, so that a — A^h^difj. Then Lemma 
[2] follows from the following argument ( |Haro , |1998[ ). In canonical coordinates, Hamilton's 
equations have the form 

9H dH 

p= ~w q= w (21) 

If we define the Hamiltonian vector field x(t) = (x q , x p ) = (q,p), then Hamilton's equations may 
be succinctly written ( Marsden fc Ratiu] 1994) in terms of a as 

(22) 



da = -dH. 



By the definition of the Lie derivative (Theorem 4.3.1 of Marsden &; Ratiu, 1994), using Cartan's 



magic formula, and the fact that differentials and pull-backs commute, we have 

^(/ 8 *«) = / s *^x( s )" = /s («x(s) d " + di x ( s )a) = d/*(i x(s) a - H). 



(23) 



Integrating both sides from s — to s — z yields the desired formula 



f*a - a = d( / /* (i x(s) a - H) dsj . 



(24) 



Taking z = 1 in this expression and recognising (19) gives 



F*a-a = d(J fZ(pq-H)dz)=d(J (pq - H) o f z dz) = dA, (25) 

which is Lemma [2j 

To show that Lemma [2] holds even in non-canonical gauge, suppose that F*a — a — dA for 
some A. Under the change of gauge A — >► A + Vx, we find a — >► a + d%, so 



F*a - a -^ F*a - a + d(F*x - x) = d(A + J^*x - X)- 



(26) 



But this is precisely dA in the new gauge (Section [2]), so if the Lemma holds in one gauge it 
holds in any gauge. 

4. Proof of Theorem Q] 

Lemma [2] allows for an elegant proof of Theorem [lj 



Necessity 

We already know that A is a necessary condition for topological equivalence, because it is an 

ideal invariant. But we can also see this from Lemma |2j Assuming F = F, it follows that 

dA = F*a -a = F*a - a = d A, (27) 



so that A and A differ by at most an overall constant. But since F\q^ = Fq^ (with the same 
winding number) we must have A — A (see Figure bj). 




Figure 3. Illustration of why F — F implies 
A — A when the winding numbers are the 
same. Both braids have the same n x A on 
dV ', although they may have different field 
lines on this boundary linking xo and F(xq) 
(as shown). But the flux through the loop 
must vanish, so A = A by Stokes' Theorem. 



Sufficiency 

To prove the converse — that A is a sufficient condition for topological equivalence — assume that 

A — A and define the mapping G — F o F _1 . Then, using Lemma bl 



G*a-a = (F- 1 )*oF*a-a (28) 

(29) 

(30) 

*<U (31) 

(32) 
(33) 
(34) 

Now we determine the possible mappings G : Di —> Di satisfying G*a = a, which |Haro| ( |1998[ 
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L )*oF*a-a 






(F~ 


l )*(a + dA)- 
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(F~ 


l )*a-a + (F- 


~ l ydA 




(F~ 


^a-iF' 1 )* 


oF*a + 


(F- 1 


(F~ 


1 )*{a-F*a) + {F- 1 )* 


dA 


-(F^ydA+iF- 


l )*dA 




0. 









2000) calls "actionmorphisms". If a is in the canonical gauge A T ^ = 0, then the possible 



mappings G that preserve a are known (Marsden &; Ratiu, 1994, Proposition 6.3.2) to take the 
form of cotangent lifts 

G(q, P )=(r(q), 1 ^\ (35) 

where T is a diffeomorphism. Since 8Dq is not a line of constant q(= -0), it follows from our 
assumption i^laoo — ^1&d that ^(?) = Q an( i hence that G g = q, G p = p. We see that G^ = tp. 
Recalling that p = A T ^h^, we note in canonical gauge that B z > implies that dp/ dp > 0. 
It follows that the coordinate transformation from (p, iji) to (p, g) has non-zero Jacobian, and 
hence that G p = p. So G = id and therefore F = F. D 

Remarks 
(i) If a is not in the canonical gauge, then G need not take the form of a cotangent lift. Indeed, 
the nature of the actionmorphisms appears to depend on the gauge of a. For suppose that 
G*a — a = for some a. Then after a gauge transformation a —> a + &x we nn d 

G*( a + d X )-(a + dx) = d(G*x-x)^0 (36) 

so that G is no longer an actionmorphism in the new gauge, 
(ii) Notice that the property of whether or not two magnetic braids are topologically equivalent 
does not depend on the choice of reference field. This choice is arbitrary; all that matters 
is that the same reference field is used to compute both A and A. On the other hand, the 
nature of differences in A between two inequivalent braids may depend on the choice of 
reference: this is an issue for future investigation. 

5. Examples: Magnetic Braids on a Cylinder 

We conclude with two simple examples on the cylinder V = {(r, (/), z) : < r < 4, — 8 < z < 8}. 
The first serves to illustrate the theory with a simple magnetic braid where A may be calculated 
analytically, and the second to illustrate an important topological interpretation of A introduced 



by Yeates & Hornig (2013) for magnetic braids in the cylinder. 



5.1. Centred Flux Ring 

Take standard cylindrical coordinates p — r and ^ = 0, so that the scale factors are h p = h z = 1, 

h^ — r. Consider the magnetic braid 

B = e z + V2re- r2 ' 2 - z2 l% (37) 
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'III 
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Figure 4. The example magnetic field from Section 5.1, showing field lines (left) and profiles 
of B ( j ) (z = 0) and A as a function of r (right). 



which consists of a localised toroidal "flux ring" in a uniform background field (Wilmot-Smith 
et al. , 2009 ) . For this field the field lines are 



f r (r, <f>; z) = r, for, </»; z) = V^e^^ 2 (erf (z/2) - erf (z /2)) + <j>, (38) 

and taking zo —> — oo 5 z —> oo m (J3TJ) , the field line mapping is 

F r (r, <j>) = r, F^r, </>) = 2 V / 2^e" r2 / 2 + <f>. (39) 

A vector potential in the appropriate gauge is 



+ V2e- r2 / 2 ~ z2 / 4 e z . 



(40) 



Since A z « on all boundaries to high accuracy, the corresponding reference is A ref = (r/2)e<p, 
giving B = e z . One can then compute the topological flux function analytically, finding 



A{r,<j>) = v^(r 2 + 2) e" 7 " 2 / 2 . 



(41) 



Notice that dA/dr — — v / 2~7rr 3 e r ' 2 and dA/d<p — 0. In this case a — (r 2 /2) d</> and we can 
verify that 



dF A 



_Fl(dF\ 

2 \dr "' ' d(p 

so that indeed F*a — a — dA. 



(**<*)( 



r, ©) " — f : -^ E dr + ^f d<fr\ = r — (-2y/tor e~ r2/2 dr + d</> 



(42) 



5.2. Displaced Flux Ring 

Suppose now that the flux ring is displaced from the origin and centred at (r = 1, 

the new vector potential is 



A = -e^> + \/2exp 



(r cos <j) — l) 2 + (r sin <j)) 2 z A 



e z . 



0). Then 



(43) 




Figure 5. Asymmetry when the flux ring is displaced from the origin (Section |5.2[ ). Panel (a) 
shows contours of A. Panels (b) and (c) show contours of the pairwise linking number c Xo ^ yo as 
a function of yo, for xo — (a/3, 55°) (in b) and xo = (\/3, —55°) (in c). The red circles show xo 
and the thick red line shows the projection of the field line f(xo',z). Solid (dashed) contours 
show positive (negative) values of c Xo ^ yo . All three panels are computed numerically. 



This is not topologically equivalent to the original example, so A must differ. One might expect 
that the A distribution would simply be translated, but Figure [5^a) shows that this is not the 
case: the profile of A within the disturbed region is no longer symmetric about the ring centre. 
This reflects the fact that A is a global quantity: it counts not just the horizontal flux of the 
ring itself, but how it is linked with the reference field B ref 
asymmetrical with respect to the position of the flux ring. 



e z . 



In (43), the reference field is 



This may be understood using the following alternative formula for A on the cylinder. |Yeates 
fc Hornig| ( 2013| ) show that A(xo) is the average pairwise linking number between the field line 
/(xq^z) and all other field lines f(yo'-,z). The linking number in this context ( |Berger , 1988) is 
the integral 



<co,2/o 



1 

2^ 



de 



xo,yo 



(*) 



dz 



dz, 



(44) 



where XQ ^ yo (z) is the orientation of the line connecting f(xo;z) and f(yo;z) in the plane at 



height z. It is proved by Yeates & Hornig (2013) that 



A(x ) = / ■ 

J D 



c X0 ,y B z (y )d 2 y () . 



(45) 



As an aside, notice that equation (45) gives a formula to compute A directly from B without 
needing to know A. 

The asymmetry in Figure [5F a) arises due to differing c Xo ^ yo for different field lines f(xo',z) 
through the flux ring. To see this, consider two field lines f(xo;z) passing through the ring, 
but starting on opposite sides (Figures [5b and [5b). The contours show c Xn _ 1ln as a function of 



^o,2/o 



7/0 for each of these xq. For yo within the ring itself, the c Xo ^ yo contours are the same in each 



case up to rotation, so these points make the same contribution to the integral (45). But for 
most yo outside the ring, c Xo ^ yo has the opposite sign in Figure J5Fb) compared to Figure 5(c), 
because of the field line direction f y (xo; z) < in Figure J5^b) while f y (xo; z) > in Figure 5(c). 
Furthermore, since c XQ ^ yo is opposite in sign on the left and right of the ring, and since there is 
a larger area to the left of the ring, the net contributions to ( [45] ) from yo outside the ring differ 
in each case. This explains the asymmetry. In Section 5.1, the area to the left and right of the 
ring was the same, hence the symmetric A profile. 
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